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MODULI SPACES OF g-CONNECTIONS AND GAP PROBABILITIES 


ALISA KNIZEL 


Abstract. We develop a g-analogue of methods introduced by Arinkin and Borodin in [T], 
|5]. Our goal is to show that the one-interval gap probability for the g-Hahn orthogonal 
polynomial ensemble can be expressed through a solution of the asymmetric g-Painleve V 
equation. The case of the q-Hahn ensemble we consider is the most general case of the 
orthogonal polynomial ensembles that have been studied. Our approach is based on the 
analysis of g-connections on with a particular singularity structure. It requires a new 
derivation of a g-difference equation of Sakai’s hierarchy m of type ^ 2 ^^. We also calculate 
its Lax pair. Following [5], we introduce the notion of the r-function of a g-connection 
and its isomonodromy transformations. We show that the gap probability function of the 
g-Hahn ensemble can be viewed as the r-function for an associated g-connection and its 
isomonodromy transformations. 


1. Introduction 

The connection between gap probabilities for orthogonal polynomial ensembles and Painleve 
equations was established in the 90’s. In continuous settings the gap probability function 
was explicitly written in terms of a specific solution of one of the six Painleve equations for 
all classical weights. It was done in [221 for the Her mite and Laguerre weight, in [ 22 ] and 
PE] for the Jacobi weight, and in m, ii for the quasi-Jacobi weight. 

The first results for the discrete case were obtained in [7j. In present paper we consider 
the g-Hahn orthogonal polynomial ensemble. We present an explit expression of the gap 
probability function for this ensemble in terms of a solution of the g-asymmetric Painleve V 
equation. The g-Hahn ensemble is the most general case of orthogonal polynomial ensembles 
associated to a a family of polynomials in the Askey scheme that have been studied in this 
context. We believe that the expressions for the gap probabilities in terms of solutions of 
Painleve equations in all known cases should follow from our result through a degeneration 
procedure. In the end of the paper we present a numerical evidence supporting this claim. 
Considering a specific limit regime we show that the distribution of the rightmost particle 
in the g-Hahn ensemble numerically converges to the Tracy-Widom distribution. 

We initially became interested in the problem of studying the g-Hahn ensemble as it 
naturally comes up in the statistical description of the tilings of a hexagon by rhombi (see 
[lOj l as well as in the representation theory (see 1211). Let us describe the former connection 
in more details. 


1.1. Tilings of a hexagon and q-Hahn ensemble. We start by recalling the definition 
|22j of g-Hahn polynomials. 
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Definition 1.1. Let q e (0,1) and N e Z>o. Let 0 < a < g ^ and t) < (3 < q ^ or a > q ^ 
and f3 > q~^. Define a weight function on X = {0,..., N} as 

{aq,q-^]q)x 


uj{x) = {a(3q)~ 




where {yi,... ,yi]q)k = {yLq)k - ■ ■ {yi]q)k, and {y]q)k = {I - y){l - yq) ■ ■ ■ {I - yq’^ is the 
q-Pochhammer symbol. 


Definition 1.2. Let us fix an integer k > 0. The q-Hahn orthogonal ensemble is a probability 
measure on the set of all k-subsets of = {j/j = : i = 0,..., N} given by 

^ ^ 

Diijxii ■ ■ ■ 1 2/xfc) = I I idJxi ~ yxj) ■ I I aj{xi), 

2=1 

where Z is the normalizing constant. 


The collection {y^^,... ,yxfi) is often referred to as fc-particle configuration. 



Figure 1. Tiling of a 4 x 3 x 2 hexagon. 

The g-Hahn orthogonal ensemble is closely related to the following tiling model. For any 
integers a, 6, c ^ 1 consider a hexagon drawn on a regular triangular lattice (see Figure [^. 
The tilings of the hexagon by rhombi are obtained by gluing two neighboring elementary 
triangles together (such rhombi are called lozenges). In Figure one can also view rhombi 
as faces of unit cubes and see the three-dimensional shape corresponding to a tiling. It is 
called a 3-D Young diagram or, equivalently, boxed plane partition. 

Consider the set of all tilings of a hexagon by rhombi, which we denote by Q{a,b,c). Let 
us denote the purple lozenges by O and introduce coordinate axes {i,j) shown in Figure 
Define a probability measure on D(a, b, c) as 

(1.1) P{T e D(a, b, c)) = where w{T) = 1~[ w{0), 

S ^^(^5) Lt 

Se% 

and w{0) = q~L°p ( is the j-coordinate of the topmost point of the purple lozenge). 

In the language of boxed plane partitions, this definition assigns to a plane partition of 
volume V (=number of 1 x 1 x 1 boxes) the probability proportional to q~^. If we send g ^ 1 
we obtain a uniform distribution on the set of all tilings. 

Consider a transformation of the hexagon by means of the following affine transformation: 
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In this way it is easy to see the bijection between the tilings and the set of non-intersecting 
paths on the plane lattice |17j . Figureillustrates this. 




Figure 2. Modihed tiling of a 3 x 3 x 3 hexagon and the corresponding family 
of non-intersecting paths. 

Let us introduce new variables T = a + b, N = a, S = c and coordinates {t,x) (see Figure 
1^. Let the parameters a, b, c be hxed. Consider the section of the modihed hexagon for 
some t. Then the coordinates of the nodes (intersection of the paths with the section) belong 
to 

= {x e Z: max{0,t + S — T) ^ x ^ min{t + N — 1,S + N — 1)}. 

Note that there are exactly N nodes for any section t. Denote the conhguration of the nodes 
by C{t). 

Theorem 1.3 ([TO]. Theorem 4.1). 

N 

Prob{C'(f) = (xi,.. .,Xn)} = const • - q~''^YY\^CN,s,T{xi), 

i<j i=l 

where Wt^N,s,T{x) is the weight function of the q-Hahn orthogonal polynomial ensemble. 

Thus, the study of the g-Hahn ensemble gives a lot of insight into the tiling model. Let 
us get back to the ensemble itself. We are interested in the following function 

^ Vigixi) ■ ■ 

which is called the one-interval gap probability function. We will call the gap probabilities 
its values for different s. Note that Dk{s) = Prob(max{xi} < s) and i/maxfxi} is called the po¬ 
sition of the rightmost particle. For the corresponding tiling model using the gap probability 
function one can describe positions of the upmost and bottommost nodes. Therefore, in the 
limit regime it provides a lot of information about the frozen boundary, a curve separating 
the so-called liquid region from the facets. 

Now we are ready to present our main result. 
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Theorem 1.4. The one-interval gap probabilities Dk{s) satisfy the following recurrence 

Dk{s)Dk{s - 2) _ a/3{rsW - g~^^+^+^)(r^+iW - 

Dk{s — 1)2 g~ 2 s(^g-s+i _ — aq){q~^+^ — b~^q-^){q-^+^ — q) ’ 


where {rs,ts) is the solution of the asymmetric q-Painleve V equation 

- b~'^q-^){ts - aq){ts - q) 


{Vs+lts + l){rsts + 1 ) = 
{rs+its+i + l)(r^+its + 1) = 

q-is+i 


af3~^q~‘^^{ts — g“^+^)2 
{q-^Vs+i + l){P~'^q~^rs+i + l)(agr^+i + l)(gr^+i + 1) 


afi{rs+iw — q '^^~^’^){rs+iw — a ^f3 ^q 2 * ^+i) ’ 

and the initial conditions can be found explicitly using Proposition 5^ and Proposition |5.6 


below. 


This result predicts the appearance of the Painleve transcendents in the limit regime. In 
Section we consider a limit regime which corresponds to a simultaneous linear growth of 
the sides of the hexagon and present a numerical evidence of the appearance of the Tracy- 
Widom distribution near the frozen boundary. We use a method suggested by Olshanski in 

| 23 |. 


1.2. Moduli spaces of g-connections. The proof of Theorem |1.4| is based on the ideas 
introduced in |2]. We start by investigating a notion of the r-function of a vector bundle on 
equipped with a g-connection. The r-function we consider is a discrete analogue of the 
continuous r-function introduced by Jimbo, Miwa, and Ueno in |18| . which proved to play a 
central role in the theory of the isomonodromy deformations. In the continuous settings the r- 
function is a holomorphic function on the universal covering space of the space of parameters 
of the connection, which vanishes when the corresponding isomonodromy deformation fails 
to exist. The g-difference Painleve equations were first studied from the point of view of 
g-isomonodromy deformations in [25]. 

Definition 1.5. Let g ^ 0,1 be a complex number. A q-connection A is a pair {C^A{z)), 
where C is a vector bundle on P^ and A{z) is a linear operator 


A{z) -. Cz -* C 


qz 


that depends on a point z e P^\{oo} in a rational way {in particular, A{z) is defined for all 
z e C outside of a finite set); here is the fiber of C over z e P^. In other words, A{z) is 
a rational map between the vector bundle C and its pullback via the automorphism P^ ^ P^ 
that sends z ^ qz. 


In speaking of g-connections, we will frequently abuse the notation and write A{z) instead 
of {C,A{z)) when there is no danger of confusion. 

Definition 1.6. ITe say that a point zq e P^ is a pole of A if A{z) is not regular at z = Zq. 
We say that zq s is a zero of A if the map A~^{z) : Cqz Cz is not regular at z = Zq. 
Note that A can have a zero and a pole at the same point. 

We proceed by associating to the g-Hahn ensemble a vector bundle on P^ equipped with 
a g-connection and constructing a sequence of its modifications. 


4 











Definition 1.7. Suppose TZ: C ^ C is a rational isomorphism between two vector bundles 
on P^. We say that C is a modification of C on a finite set S' c P^ if TZ{z) and TZ~^{z) are 
regular outside S. We call C an upper modification of C iflZ is regular {resp. C is the lower 
modification of C). 

A q-connection {C,A{z)) induces a q-connection {C,A{z)), which we also call a modifica¬ 
tion. 


The next step is computation of “the second logarithmic derivative” of the r-function for 
the sequence of modifications and proving that it coincides with “the second logarithmic 
derivative” of the gap probability function. 

In order to follow our plan, we study a special class of ^-connections in the spirit of 
More precisely, we construct moduli spaces of g-connections with a particular singularity 
structure. The choice of the singularity structure we consider is dictated by the singularity 
structure of the g-connection associated to the g-Hahn ensemble. 

Remark 1.8. After we choose a trivialization of C restricted to = F — {oo}, the operator 
A(z) can be written in coordinates as a matrix-valued function A{z). ITe will call it a matrix 
of the q-connection. For two different trivializations the corresponding matrices differ by a 
q-gauge transformation 

(1.3) A{z) = R{qz)A{z)R~^{z). 

Thus, the moduli space of q-connections can be identified with the moduli space of their 
matrices modulo q-gauge transformations. 


Let us consider g-connections {A{z), 0@ 0{—l)). We impose a number of restrictions on 
them. After choosing a trivialization the matrices of the g-connections are 2x2 matrices 
with polynomial entries of the following form; 


A(;.) 


Oil ai2 
^21 ^22 


A(0) 


w 0 
0 w 


deg{aii) A, 3, deg{ai2) < 4, deg{a2i) ^ 2, deg{a22) ^ 3 
det(A( 2 ;)) = uv{z — ai){z — a 2 ){z — a 3 )(z — a 4 ){z — a 5 )(z — Qq), 
where ai,.. .Oq, u, v, w are complex parameters. Denote 

Q( \ r 1 0 1 

S{z) = ^ ^ . 

z _ 

We also require that 

det{S{qz)~^A{z)S{z)) = quvz^ -f 0{zA and tT{S{qz)~^A{z)S{z)) = {u -\- qv)z^ -f 0{zA- 


Note that S is essentially a basis of O @ 0{—l) in the neighborhood of {oo} e P^. We say 
that such g-connections are of type A = (oi,..., og; n, gn; w, w; 3). 


We consider g-connections modulo g-gauge transformations of the form (1.3). We can write 
R in coordinates as 


R(z) 


^11 ri2 

0 r22 


deg(rn) ^ 1, deg(ri 2 ) < 2, deg{r 22 ) ^ 1- 
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Lemma 1.9. Under certain nondegeneracy conditions on the parameters ai,... ,aQ,u,v,w of 
a q-connection A, there exits its unique modification A of type {qai, qa 2 ... ,ae;u,qv; qw, qw;3). 


Let us assume that the parameters oi,..., ag, u, u, w are generic. We will discuss what it 
means for the parameters to be generic in Section We proceed by showing that the 
moduli space Mx of g-connections of type A = {ai,... aQ;u,qv;w,w]3) modulo g-gauge 
transformations is two-dimensional and its smallest smooth compactification can be identified 
with blown-up at nine points; more precisely, it is a Sakai surface of type ^4^^^ Let us 
state our next result. 


Theorem 1.10 (g-Py). Consider the second root of the matrix element 021 {note that 0 is 
always a root 0 / 021 ) os the first coordinate on Mx, denote it by t. Let the second coordinate 
be 

w{t — afi){t — afi){t — a^){t — Qq) 1 
011(^)03040506^ t 

Consider the modification of C to C from Lemma 


1.9. It shifts 


oi ^ goi, 02 go 2 , w qw. 

Then this modification defines a regular morphism qPV between two moduli spaces Mx and 
M^, the moduli space of q-connections of type 

A = (goi, go 2 ,..., Oe; u, go; qw, qw, 3). 


Moreover, the coordinates {t, r) on the moduli space M^ are related to {t, r) by the asymmetric 
q-Painleve V equation 


(1.4) 


{rt + l){rt + 1) 
{ri + l){ri + 1) 


uvafaKras + l)(ro4 -I- l)(ro5 -I- l)(ro6 -I- 1) 
{rw — vaia 2 ){qrw — UO 1 O 2 ) 

QiQ2(t - 03)(t- 04)(t- 05)(t- Qe) 
03040506(gf — Oi)(gf — O 2 ) 


In the proof we also explicitly compute a Lax pair for this equation in a new fashion. 

The asymmetric g-Painleve V equation was first found by B. Grammaticos, A. Ramani 


and their co-authors [Uj. The form of the equation 1.4 we use differs from the standard one 
by a change of notation. 


In [26] Sakai showed that the asymmetric g-Painleve V equation 
appears as a particular case of the four-dimensional g-Garnier system {N = 2) and first 
expressed it in a Lax formalism. The computations done by Sakai are not applicable in our 
situation since the modifications of the parameters we aim to study are different from the 
ones he considered, see Section for details. The other two approaches to obtaining Lax 
pairs were presented in [23], |31j . 

We believe that g-orthogonal ensembles, for which the similar results were obtained in [7], 
also fit in the framework of g-connections. 


The paper is organized as follows. 

In Section we discuss the properties of g-connections and show that the geometric ap¬ 
proach to isomonodromy deformation of g-connections implies that the transformations lift 
to isomorphisms between surfaces. We also define modifications of a special class of g- 
connections which will lead to the asymmetric g-difference Painleve V. 

Then in Section]^ we give the proof of this result. 
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Section is devoted to the discussion of the g-analogue of the constructions introduced in 
[2]. The g-generalization is mostly straightforward, however, the computations require some 
work. We conclude this section by providing the recipe for computing the r-function of the 
isomonodromy transformations leading to the asymmetric g-Painleve V equation. 

In Section we investigate the g-Hahn orthogonal polynomial ensemble and give a proof 
of Theorem 11.41 

In Section we present a numerical evidence of the appearance of the Tracy-Widom 
distribution near the frozen boundary of the tiling model of a hexagon by rhombi with the 
weight proportional to 

Acknowledgements. The author is very grateful to Alexei Borodin for suggesting the 
project and for many helpful discussions. 

2. Moduli spaces of g-coNNECTiONS 


2.1. Preliminaries. 

Definition 2.1. Let C he a rank m vector bundle on and A{z) be a q-connection on C. 
Suppose A(z) satisfies the following conditions: 

(i) The only zeroes and poles of A{z) are as follows: a pole of order n at infinity and 
simple zeroes at k distinct points Oi,..., oa, e (he., A(z) is regular and det(A( 2 ;)) 
has zero of order one at Oj). 

(a) On the formal neighborhood o/oo e P^, there exists a trivialization TZ{z) \ LA ^ C such 
that the matrix A with respect to IZ satisfies 

n{qz)A{z)n-\z) = Anz'^ + oiz'^-^), 
where An is a semi-simple matrix with eigenvalues pi,..., pm- 
(Hi) On the neighborhood of zero A{z)\z=o = A^ is semi-simple with eigenvalues 6 *i,... Am- 

We call q-connection A{z) or, more precisely, the pair {C,A{z)) a q-connection of type 

A (®1) ■ ■ ■ 1 Ply ■ ■ ■ y Pmy ^1) • • • ) ^my P') • 

Remark 2.2. ITe will also consider q-connections that have simple poles besides simple 
zeroes. The operation “multiplication by scalar” turns a pole into a zero and vice versa. 
More precisely, let f{z) ^ 0 be a rational function on P^, and let A{z) be a q-connection on 
C. Then the product f{z)A{z) is again a q-connection on C. For any q-connection A{z), we 
can choose f{z) such that the only pole of the product f{z)A{z) is at infinity. In this case 
we will slightly modify the notation for the type of the connection. We will write A is of 
type (ai,..., afc; 6 i,..., bp pi,..., pm] 6 * 1 ,..., dm] n), where (oi,..., a^) and ( 6 i ,... ,hi) is the 
collection of zeroes and poles respectively. 

Definition 2.3. Let A{z) be a q-connection on C of type (oi,..., 0 ^; pi,..., Pm] Oiy ■ ■ ■ y Om] n). 
The q-degree of A{z) is the following quantity 

1 —rk 

degg(A(z)) = '"nm 

1 h=i 

Lemma 2.4. Let A(z) be a q-connection on a rank m vector bundle C of type 

(fll, . . . , OA: 1 P\ y ■■■ y Ply b^ly ■■■ y j 
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with invertible Aq and An- Then the following holds 

mn = k — I and = degg(^( 2 :)). 

Proof. Determinant of the morphism is a map det(^( 2 ;)) : /\f^ Cz h^qz- In 

this way, we can dehne a g-connection det(^) on a line bundle det(£). Note that det(^) is 
of type (ai,..., a^,; Y[T=i Pk OLi Oi]k-l). □ 

Corollary 2.5. Suppose A{z) is a q-connection on rank 2 vector bundle C of type 

• • •, U/;,, 6]^,..., 6;, p2; ^5 n) 

with invertible Aq and An- Suppose that for any I a {1,... ,k} we have A for any 

iel 

k E Tj and j = 1,2. Then A{z) is irreducible, i.e., there is no rank one subbundle i cz C such 
that A{iz) zz iqz for all z. 


Proof. Suppose that £ e £ is an invariant subbundle of rank 1. Then A induces a g-connection 
Ai on 1. All its zeroes belong to {ui,..., a^} and its type is either (oi,..., a^; pi; Q\ n) or 
(oi,..., Ofc; P 2 ; d-,n). Now Lemma 2.4 leads to a contradiction. □ 


Lemma 2.6. Suppose that A is an irreducible q-connection on a rank 2 vector bundle C of 
the type A = (oi, ..., a*,; pi, P2; 6*1, 6*2; u) with invertible Aq and An- If C ^ 0{ni) ©0(^2), 
then \ni — n 2 \ ^ n. Moreover, if qpi = P 2 the inequality is strict. 


Proof. Without loss of generality we can assume that ni ^ 77 . 2 . Let £ c £ be a rank 1 
subbundle of degree ni. Since (£, Al) is irreducible, i is not invariant, thus, a rational map 
fj-. i. C ^ s*C s*CII is not identically zero (recall that s\ z e qz e P^). 

Notice that can have at most an order n pole at 00 (and no other poles). Therefore, 
Hi = deg(f') < n + deg(£/£) = n + 712 - The second statement follows from the fact that if 
gpi = p 2 then has at most an order n — 1 pole, because the coefficient of z"^ in is an 
off-diagonal element of a scalar matrix. □ 


From now on we will work with g-connections on rank 2 vector bundle. Denote by Mx the 
moduli space of g-connections of type A = (oi,..., a^,; &i,..., bp pi, p 2 ] 0, 9] 2). From Lemma 


2.4 we see that Mx is empty unless 


(2.1) k = 2n, 

(2.2) deg^(A) = g^, where k is an integer . 
Let us also consider the following assumptions on A: 


(2.3) 

(2.4) 


«n 


iel' 
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A g^ for any / c {1,..., fc}, any k e Z and j = 1, 2; 
Pi, P 2 A 0 9i = 92 A 0. 


Assumption (2.3) implies irreducibility of the g-connection and can be used to prove that 


the moduli space Mx is a smooth variety. In the present paper we do not work with the 
categorical dehnition of the moduli space since it is not necessary for our purposes. However, 
we think that one can prove that the moduli space we consider is the coarse moduli space 
following the arguments in [3|. 

In what follows we will need more facts about the g-connections. 






Definition 2.7. Let be the set of all q-connections satisfying 2.1-2. f 


Definition 2.8. Let C he a rank two vector bundle on a: e P^ and I Lx a one¬ 
dimensional subspace. Denote by L the sheaf of sections of C. An elementary lower modi¬ 
fication of C at X with respect to i is a rank two bundle C whose sheaf of sections is 
L = {s e L\s{x) E i}. 

We can dually define elementary upper modifications. 


Let us denote the set of poles and zeros of A by Sing(.4,). 

Lemma 2.9. (Lemma 1.7 P]) Suppose that q~^x is not a singular point of A, while x is a 
singular point of A. Then there exists a unique modification such that 

(i) Sing(Al^®^) = (Sing(A) x) u q~^x] 

(%%) A is the unique modification of A^^^ at x with no singularity at q ^x. 

The dual statement also holds. 

Lemma 2.10. Suppose {C,A) e Mx for X = X = (oi,Ofc; pi, p 2 ; 6^, 6*; n) e A^. Let {C,A) 
he an elementary upper modification of C at x eF^. Then the only cases when (£, A) belongs 
to for some A e A^ are as follows: 

(i) If X = cc, then X = (oi, • • • ,ak;q~^pi, p2;0,9;n) or X = (oi, ...,ak]pi,q~^p2]9,9'A)- 
(a) If X = Oi is a zero of A and q~^x ¥= aj is not, 

then X = (oi,--- ,q~^ai,--- ,ak',pi,qpi,9,9;n). In either case, the elementary modifi¬ 
cations define an isomorphism Mx —» M^. 

2.2. Geometric description of Mx- We will focus on g-connections of type 

A = (oi, ...,ak] u, qv; w, w; 3) e A'’ 

on C ^ O This choice is dictated by appearance of this type of connections in our 

analysis of the g-Hahn ensemble. 

Let us choose an isomorphism iS; £ ^ O Then A induces a g-connection of type 

A on (9 0 0{—l). Such a g-connection is represented by a matrix of the following form 

^ _ Oil ai2 

[ 0.21 022 

This representation is not unique: S can be composed with an automorphism of the 
bundle. Such an automorphism can be written as a matrix 


an, a22 e r(P\ 0(3)), e T{F\ 0(4)), a2i e r(P\ 0(2)). 


(2.5) R= '■(5' , rn,r22eC-{0}, m e r(P\ 0(1)). 

If we choose another trivialization S then A is replaced by the g-gauge transformation 


R{qz)A{z)R{z) 


-1 


Lemma 2.11. Let A be a q-connection on 0©0(—1) and let its matrix be of the form (4.1) 
Denote 

" 1 0 


S{z) = 


0 i 
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We claim that A is of type X = {ai, ...,ak;u,qv; w, w; n) if and only if A satisfies the following 
conditions 

detyl(2:) = uv{z — ai){z — a2)(z — a3){z — ai){z — a5){z — Og), 
det(S'(g^)“^y 4 ( 2 ;)S'(^)) = quvz^ + 0 {z^) and tT{S{qz)~^A{z)S{z)) = {u + qv)z^ + 0 {z‘^), 

y 4 ( 0 ) has an eigenvalue w of degree 2 . 


We can now think of Mx as a quotient of all matrices that satisfy Lemma 2.11 modulo 
g-gauge transformations (2.5). 

Remark 2 . 12 . In | 26 | Sakai considered the matrices of the form 

A(^z^ = Aq + Aiz + A2 z‘^ + A'jZ^, 


where 


^3 


Ki 0 

0 K2 


det2l(2;) = k,ik,2{z — ai){z — a2){z — a3){z — a4){z — a5){z — Qq). 


He also imposed the condition qni = K2- Sakai considered the modifications that shift 

ai qai, 02 qa2 and 61 -> qOi, 62 q62, 

where 6*1, 62 are the eigenvalues of Aq. 

Theorem 2 . 13 . The smallest compactification of Mx is isomorphic to the surface of type 
^2^^ in Sakai’s classification. 


Proof. Denote by f e the second root of 021 (note that zero is always a root of 021 ). From 
the irreducibility assumption it follows that this coordinate is correctly defined. The direct 
computation shows that t is invariant under g-gauge transformations. 

It is convenient to reduce the matrices of a g-connection A of type A to some normal 
form. Solving a system of linear equations for the coefficients of R we can reduce A{z) to 
the following form (we assume f e P^ — { 00 }) 


sz + w h^z^ + + h2z‘^ + hiz 

z{z — t) k^z^ + k2z‘^ + kiz + w 


Then s can be considered as the second coordinate in this chart. 

Since A is of type A, we can find all the coefficients of its matrix in terms of s and t. We 
get 

( 1 ) hi{ts + w) = Fh^{t,s), 

( 2 ) h2 = Fh^{t,s), 

(3) /i3 = Fh^{t,s), 

(4) ki = Fk,{t,s), 

(5) k 2 {ts + w) = Fk^it^s), 


where F_ is a polynomial. Moreover, 


^ , ui. 




and 


^ / w. 


F 


Thus, in order to resolve the singularities we need to blow-up in six points (a*, 


10 










Let us work with t in the neighborhood of oo. Let us denote ^ by x. In this chart we can 
reduce the matrix of the g-connection to another normal form 

_x yz^ + w + h'iz'^ + h2Z^ + hiz 

^ ' z{l — xz) k^z^ + k2z‘^ + kiz + w 

Let y be the second coordinate. Again, since A is of type A, we can find all the coefficients 
of its matrix in terms of x and y. We get 

( 1 ) hix{x^w + y) = Gh^{x, y), 

( 2 ) h 2 {x^w + y) = Gh^{x, y), 

( 3 ) h^{x^w + y) = Gh^ix^y), 

( 4 ) ki{x^w + y) = GkAx,y), 

( 5 ) k 2 x{x^w + y) = Gk 2 {x,y), 

where is a polynomial. 

Moreover, 

Gh^O^y) = wq-^{qv - y){u - y), Gk^{0,y) = q~^{qv - y){u - y). 

Notice that in the neighborhood a; = 0 (f = oo) the coordinate y can not be zero due to our 
assumptions about the behavior of the g-connection at oo. Therefore, in order to resolve the 
singularities we need to blow-up in two points ( 0 ,gn) and ( 0 ,m). 

We covered Mx by two charts x — {oo} and x P^ — {0} blown-up in a number of 
points. We can glue them together along their intersection g 7 ^ 0 , oo. The transition map is 
y ^ s = yt^. 

Let us recall that we would like to introduce the following coordinate on Mx 

w{t — az){t — ai){t — a^){t — ao) 1 

We see that the rational map given by this coordinate r: Mx —» P^ is actually regular, since 
we have already resolved its singularities. 

Summing up, we have realized Mx as an open subscheme of the blow-up of P^ x P^ in 
eight points; more precisely, we can check by direct computation that its compactification is 
isomorphic to □ 


3. Proof of Theorem inni 


3.1. q-Difference Pv In this section we discuss the proof of Theorem |1.10[ It is based on 


calculations which are hard but straightforward. The trickiest part is the introduction of the 
coordinate r. 

Proof. Take (A, A) e Mg and set (A, A) = qPV(A, A). Note that it suffices to check the 
formulas ( |1.4 ) on a dense subset of Mg] we can therefore assume that t((A,A)) ¥= 0, 00 . 
Then there exists an isomorphism S : O ^ C such that the matrix of A relative to 

S is of the form 


(3.1) 


A{z} 


sz + w hiZ* + hsz^ -1- h 2 z‘^ + hiz 
z{z — t) k^z^ + k 2 z‘^ -f kiz + w 
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Let us introduce the coordinate p such that s = 2^. An analogue of this coordinate 
was used in |25], jT], jT9], |26|. The specihcation of the behavior at inhnity implies that 
= qu + w and = —uv. Moreover, 

6 

det(A( 2 :)) = uv{z — ai){z — a2){z — a3){z — a4){z — a5){z — qq), uv a* = w'^. 

This information allows us to express h^, ^ 2 , hi and A; 2 , as rational functions in t and p. 
By the dehnition of qPV, the matrix A is the g-gauge transformation of A: 

A{z) = R{qz)A{z)R~^{z), 

where R is of the form 


such that 


R(z) 


diZ + Og 1^2^^ T T /dg 
7o hiz + So 


det(-R( 2 ;)) = c{z — ai){z — 02 ). 


Also, we know that A has no singularities at Oi and 02 - These restrictions yield polynomial 
equations on the coefficients of the gauge matrix. The resulting system determines R up to 
a multiplicative constant and we obtain expressions for t and p in terms of f, p. The hnal 
step is the introduction of the coordinate r. 


We knew form the geometric considerations (Theorem 2.13) that we should obtain the 
asymmetric q-Painleve V. Thus, we expected the existence of the coordinates such that the 
formulas can be written in a nice multiplicative form. We investigated properties of the 
analogue of r in the case of d-connections and then came up with the formula for r in the 
case of g-connections. First, we were able to express t as a rational function in t and p. The 
next step was to factor the denominator of this expression and compute the residues. After 
some work we obtained the right candidate for the coordinate r. 


□ 


Remark 3.1. Note that A{z) = R{qz)A{z)R ^{z) is the compatibility condition for the 
following system of equations 

Y{qx) = A{x)Y{x )— linear q-difference equation, 

Y{x) = R{x)Y{x) — the deformation equation. 

Thus, we have expressed the asymmetric q-Painleve V equation in a Lax pair formalism. 


4. Computing the ratios of the t-function 

4.1. r-function. Following [2], we introduce the notion of the r-function for a vector bundle. 

Let det(C) denote the top exterior power of a hnite dimensional vector space V and 
denote its dual. Let us also dehne a one-dimensional vector space 

detRF(£) = det(R°(P\ £)) (x) (det(H^(P\ C)))-\ 
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Definition 4.1. Suppose that a rank m vector bundle C is of degree —m. 
detRr(£)-i by 


I 1, £ ^ 

[ 0, otherwise. 


We define r(£) e 


By itself r(£) contains almost no information, let us show how we can interpret the ratios 
of the r-function in some specific situations. 

Example 4.2. Let us consider a vector bundle C with slope —1 and its upper modification 
C such that H^(P^,£) vanishes. Then the long exact cohomology seguence for 

0 ^ ^ ^ £/£ ^ 0 

provides an identification detRr(£) = det(il°(P\£) (8) (det(H^(P^, £/£)))“^. If we now 
consider the morphism 

n: H°(P\£) ^H°(P\£/£), 
then t{C) can be identified with det(n). 

Remark 4.3. Since the slope of C is —1, we see that H°(P^,£) = 0 if and only if C = 

Example 4.4. Let £2 be a modification of C\. We assume that £1 is isomorphic to — 
and £2 is of slope —1. One can always find an upper modification £ of £1 that is also an 
upper modification of £ 2 . Then the ratio 

e detRr(£i) (g) detRr(£2)"^ 

corresponds to the determinant of the composition 

H°(P\£/£i) ^ H°(P\£) H°(P\£/£2). 


We will see below that if the vector bundles are equipped with g-connections then the 
second ratio ("second logarithmic derivative") of the r-function makes sense as a number. 


4.2. General q ^connections. 

Remark 4.5. We will be using q~^-connections in the rest of the paper to make all the 
computations parallel to the ones done in j2]. 

Let £ = C>(—1)”^ be a vector bundle and .4, be a g~^-connection on £. Suppose that A 
has singularities at n distinct points Oi, ... ,an and no singularities at for fc e Z — {0}, 
i = 1,... ,n. We impose no restrictions on behavior of A elsewhere. 

Consider 


fl: = = (g “Wi,...,g = (mi, ..., e Z” - {0}} c 

Let au e fl. We claim that there exists a unique modification £„ of £ such that 

Sing(A) = Sing(4,)\{ai,..., 0 ^} u {g"“"ai,..., g"“"a„}. 

This can be done by a sequence of elementary upper and lower modifications. Let us set 

: = detRr(£„)-i (x) detRr(£„_eJ. 
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Proposition 4.6. For any u, v e fl there exists a canonical isomorphism Su'^ ^ Sv\ 

Proof. The construction and the proof tautologically repeats |2j (Proposition 2.1). □ 


Let us now define the ratios of r. Set 

n 

Z” = {f = (fi,..., I 2 KiSi = 0}, 

where Hi are the degrees of the singularities a* of A. 

For any a„ e hi such that £„ is of slope -1 we see that Cu+s is also of slope —1 for 
s = (si,..., Sn) e Zq. Once t is fixed, “the second logarithmic derivative” of the r-function 

makes sense as a number. Here we used Proposition |4.6| to be able to drop the subscript of 


'^(^u+t+s) 

t{Cu + t) t{Cu+s) 


4.3. Computations in coordinates for connections with simple zeros. Now let us 

make the construction above more explicit. All the computations are parallel to the ones in 

m- 

We will present the result for one special case which we will use later. 

Let us fix an isomorphism C ~ 0{—\) © 0{—l). Then the g“^-connection A is given by 
its matrix A{z). Suppose A has singularities at Oi,..., and no singularities at q^Oi for any 
k E — {0}. The singularities at Oi,..., otv are simple zeroes and poles. Let Oi be a simple 
zero and 02 a simple pole. 

For generic A{z), there exists a unique rational matrix R{z) with the following properties 

• A{z) = R{q~^z)A{z)R~^{z) has the same singularity structure as A{z) with singu¬ 
larities at hi = gai, 0,2 = qa 2 , = a^,... An = CLn] 

• R{ao) = /; 

• Both R{z) and R~^{z) have singularities only at Oi and 02 - 

It is not hard to construct such a matrix R{z). Let us choose a basis w in the kernel A(ai) 
and a basis w' in the image of ReSa 2 A^(^), by we denote the matrix transpose. Then 

R{z) = / H-—, R^^{z) = I - - —, det(i?(z)) = - 

z — 02 z — ai z — 02 

{02 - Oi)w ■ 

-no - -7-. 

{w, w ) 

The choice of w and w' for all singular points of A{z) determines bases in the corresponding 
spaces for all deformations. More precisely, if A{z) is the next deformation such that 

(hi, 02 ,..., On) ^ (ghi, gh2, • • •, a^), 
w = R{qoi)A~^{qoi)w, w' = w''^ A{qo2)R~^{qo2). 

Now the second ratio of the r- function is equal to 

{w, w'} 

{w, w'} 
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and does not depend on the choices we made. 

Theorem 4.7. Consider q~^-connections on C = 0{—l) 0 0{—l) of type 

A = (ai, Os, Os; 02 , 04 , uq ; u , q~^v; w, w; 3). 

Then their matrices can be written in the following form 

1 


^ (z - a 2 )(z - a 4 )(z - ae) 
Also, the following holds 

det(y4(2;)) = uv 


uz^ + 92 ^'^ + giz + w 

l 2 z‘^ + l\Z 


k2Z^ + kiz 

vz^ + h 2 z‘^ + hiz + w 


{z - ai){z - a 3 ){z - as) 
{z - 02 ) (z - 04 ) (z - ae) ‘ 


We consider such matrices modulo q-gauge transformations of the form (|2.5|), where 

- f f z 


rii = const, deg{ri2) ^ 1 , ^"22 = const. 

Consider the second root of the matrix element 021 {note that 0 is always a root of 021 ) as 
the first coordinate on Mx, denote it t. Let the second coordinate be 

w{t — a 3 ){t — a 4 ){t — a 5 ){t — ae) 1 
011 (^) 03 ( 340506 ^ t 

Consider the modification of C to L that shifts 

01^01= goi, 02 ^ 02 = qa2, w ^ w = qw. 

Then the coordinates {t, f) on the moduli space are related to {t, r) by 
if we consider the modification of C that shifts 

oi —> goi, 02 qa2, w qw 

the “second logarithmic derivative” of the r-function makes sense as a number and the fol¬ 
lowing holds 


(1.4). Moreover, 


(4,2) 


2 {qrw — vaia2){qrw — ua\a2){qt — ai){t — 02) 

13 T = -- = -. 

(t(£))^ («i - '?a3)(<3i - go5)(o2 - 904)(02 - qae)aia2 


Proof. The first part of the theorem follows from Theorem 1.10 Note that there is a natnral 


modification 0 {—l) 0 0 (— 1 ) —>(9 0 0 {—l) snch that the corresponding parameters of the 
g-connections change accordingly. 

Next, to compnte the “second logarithmic derivative” of the r-fnnction we jnst have to 
follow the algorithm described above. □ 
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5. Gap probabilities and the t-function 


5.1. The q-Hahn ensemble. Let q e (0,1) and N e Z>o. Let 0 < a < and 0 < /d < 
or a > q~^ and /d > q~^. Fix a positive integer k. Recall that the weight function for the 
g-Hahn ensemble is given by 


u{x) = {a/3q) 


-X ;?) 


(g,/3-iiV-i;g)^’ 

where (i/i,.. .,yi]q)k = (l/i;g) • • • {.yi]q)k and {y,q)k = {I - y){l - yq) ■ ■ ■ {I - yq^~^) is the 
g-Pochhammer symbol. 

We consider a probability distribution on the set of all /c-subsets (often called particle 
conhgurations) of 

X = {i/, = g-G* = 0,...,iV} 

given by 

. k 




-xi _ q-xjy 






i=l 


where Z is a normalizing constant equal to 2 • • • > yxk)- 

Vxi 

Definition 5.1. Given a particle configuration {yxi, ■ ■ ■ ,yxk) we will call i/maxfxi} the right¬ 
most particle and ymin{xi} the leftmost particle respectively. 


We are interested in the gap probabilities 

-Dfc(-s) = ^ 7^{yxi, ■ ■ ■, yxfi■ 

Xl,...,Xk<S 

We will need the following fact from the general theory of the orthogonal polynomial 
ensembles. Consider a kernel /C on X x X dehned by 


yj) 


yi - yj 

u{i){u{yi)v{yi) - u{yi)v{yi)), 


i 


i = j. 


Here m(C) = Pk{C) and v{C,) = {Pk-i, Pk-i)^^^ • -Pfc-i(C)) where Pk is the monic orthogonal 
polynomial of degree k associated to the weight function and 


{f{C)9{C))w- = J] fiyx)9{yx)uj{x) 

VxSX 


is the corresponding inner product on the space of all polynomials. 

Let us denote {ys,ys+i, ■ ■ ■ ,yN} by 91^. It is well known that the gap probability can be 
expressed as a Fredholm determinant 


Dfc(s) = det(l -/Cs), s E Z^k, s ^ N, 


where ICg is the restriction of K, to 91s x 91*. 
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5.2. Geometric construction. We will closely follow the exposition in [2] (Section 5) adapt¬ 
ing it for our specific needs. Let us consider on a vector bundle 

£0 = 0{k-l)®0{-k-l). 

For any subset c X, define the following modihcation: 

• £ 9 jt and C coincide on P^\fOT; 

• Near any y e VJH sections of are rational sections s = (si, S 2 )* e £, such that S 2 is 
regular at y and Si has at most a hrst order pole at y with Resy(si) = uj{y)s 2 {y). 

Lemma 5.2 ([2], Proposition 5.2). Under the above assumptions Cx — 0{—l) 0(P(—1). 


Let be a modihcation of on — 971 whose sections near x e X — VJt are of the form 

s = (si, S 2 ) e Cxn, 

where Si is regular at x, S 2 has at most a hrst order pole at x and Res^^s '^ 10 ^= . Note 

that is also an upper modihcation of Cx- 

For every point x e X — iM consider two functionals on sections of 

fxis) = {w'x + w"x{z -x))s\:,=x, 


gx{s) = Res^s • w^, 

. Note that the sections of Cx are exactly sections on 

on which fx (resp. px) vanish for all x e X — 91T. In this way, assuming that |X| is hnite, 
one gets the identihcations 

fx-m = if.Ux-m-. ^\^\CZ/Cx) - 


where w' = 


u:{x) 

0 


and w" = 


9x-m — {.9x)xex-m.'- H°(P^,£^/£an) ^ 

This induces an isomorphism 

detRF(£x) ® detRF(£a„)-i = det(H°(P\ £“^/£x)) ® det(H°(P\ = C. 

Thus, one can interpret the ratio as the determinant of the composition 

Cl^l-I®^l ^ H°(P\£^/£-e) ^ H°(P\£^) ^ H°(P\£^/£ot) 


Theorem 5.3 ([2], Theorem 5.3). For any 9Jl c X, the following holds 

(5.1) ^ det(l — K,\p{x-m))- 

'TyCx) 

Proposition 5.4. Consider £0 = 0{k — 1) © 0{—k — 1) and TUs = {1, • • •, ^ X, 

where s ^ k. Then £®t — 


Proof. From the construction it follows that deg(£fOT) = —2, thus, it is enough to show that 
£ot has no global sections. 

The section of C^a is of the form s = (si, S 2 ), such that Si is of degree at most s — 1 and 
S 2 is given by 
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( 5 . 2 ) 


s—1 

i=0 


w{q *)si(g *) 
2 ; — 


and its order of zero at = 00 is at least k + 1. 

Equivalently, it means that for any polynomial p{z) of degree at most k — 1, 

s-l 

lim p{z)s 2 {z) = y.piq~')w{q~')si{q~^) = 0 , 

•?:—».rr) < * 


2 = 0 


which is impossible. 


□ 


5.3. q ^-Connection associated to g-Hahn ensemble. Consider the following g ^-connection 
on £0 


<?^(<? g 

u]{z) 

0 1 


24o(^) = 

Then the corresponding g"^-connection on £ 531 ^ is given by 


{z-aq)-(z-q p, 

ctl3{z—q)-{z—b~^q~^^) 

0 1 


As{z) = Ms{q ^z)Ao{z)M,\z), 

where M.k{z) is the unique solution of the following normalized Riemann-Hilbert problem; 
• It is analytic in C — {1,... 

0 u:{i) 


Res M-kiz) = lim M.i 


z=q 


z^q 


0 0 


Mk{z) 


z-^ 0 


0 z^ 
In [7] it was shown that 


= I + O(-), as 2; —> 00. 


Msiz) = 


Pkiz) 


y PkiAAA 

^ Z—U 

uG{l,...,q-^+'^} 

cP,-,{z) s ayya 


where Pk is the monic orthogonal polynomial of degree k with respect to the g-Hahn weight 
on {1,..., g"*+i} and c = (Pfc-i, Pk-i)Z^- 

Let us compute Ak{z). By substitution one can check that 


Mk{z) = 


k—1 

n (^ - (i~") 

2 = 0 

k—1 k—1 k—1 

n (^ - < 1 ~") 2 Y\{z-q 

2 = 0 2 = 0 2 = 0 


— 2\ —1 


where p* = ^ U 


O^j^k—l 


q-‘-q- 
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Thus, 


Ak(z) = 


{z — q~’^+^){z — /3~^q~^){z — q) 

3^{z - q){z - q-^){z - aq) 0 

021 ( 2 ;) 


where 


- ^—-{z 


k-l 


k-1 


-aq){z-q - Q) Yj q-i^ 1 q-m 


z — q~ 


i=0 ^ i=0 

Let us denote Oi = 02 = 03 = 04 = (3~^q~^; 05 = aq; Qq = q. 

Summarizing, we get the following statement. 


Proposition 5.5. The matrix of the q ^-connection Ak{z) is of the following form 


Ak{z) = 

with 


{z - a 2 ){z - 04 )- og) 

det(Afc(2;)) = 


'^z^ + g 2 Z^ + giz + w 


0 


f 2 z‘^ + fiz q^z^ + h 2 Z^ + hiz + w 

{z - ai){z - a^){z - Os) 


0/3(2; — a2)(2; — a 4 ){z — qq) 

Introducing the usual coordinates t and r, we can directly compute 

• w = j3~^q~^~^'^‘^; 

• Tfc = -04 

• tk = where 

*^-1 q-k+i 

fi = Yj ~ + ® 3®5 + ® 3®6 + 05 ^ 3 ) 


m=0 


k-l 


al3 

— Q^{.Q — (^1 ^2 ^4)0' ^ + ( i \( X2 " 1 “ ( X \( i ^ + (^2^4))! 


Q — k-\~l 

f 2 = Y - 03 - 05 - Oe) - q^{q~'^ - Oi - 02 - 04)). 


m=0 


afl 


Notice that the isomonodromy transformation .4,(s) ^ .4,(s + 1) shifts 
Oi —> Oi = goi, 02 ^ 0,2 = qo 2 , w ^ w = qw. 


The only thing left in order to apply Theorem 4.7 to the series of modifications ^^( 2 ;) and 
prove Theorem 1.2 is to compute the initial conditions, i.e., Dk{k), Dk{k + 1). 

Proposition 5.6. ([7j, Proposition 6.6) Let'X. c M 6 e o discrete set, let {P„( 2 ;)} be the family 
of orthogonal polynomials corresponding to a strictly positive weight function a; : j£ ^ M. 
Then 


k-l 


(5.3) 


Dtik)--- [~[ (it,- ir,)^ ■ ]/[w(;), 


O^i^j^kl 


1=0 
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k-1 


(5.4) Dk{k + 1) = oj{k) ■ Dk{k) ■ ]^(7rfc - Wif, 

where qk is given by 


1=0 




Qk = 


{pk + ^ 


Pn 




)• 


7n=0 


k-l 


Remark 5.7. Using the standard arguments one can show that Z is equal to ]))[ (Pj, Pj)^ 

i=0 


We are finally ready to prove Theorem 1.2 


Proof. What we need to show is that the second logarithmic derivative of the r-function 
for the series of isomonodromy transformations As{z) coincides with the second logarithmic 


derivative of Dk{s). If so, the theorem will follow from Theorem 4.7 


The ^-connection Aq(z) was chosen in such a way that we could apply Theorem 5.3 to 
obtain the result, see j2]. □ 


6. Numerical computations 

The goal of this section is to show the numerical evidence supporting the conjecture, stat¬ 
ing that the distribution of the rightmost (resp. leftmost) particle in the g-Hahn ensemble 
converges to the Tracy-Widom distribution after proper centering and scaling. In the cor¬ 
responding tiling model this describes fluctuations around the frozen boundary. We do not 
have the rigorous proof of this conjecture, but we expect that it can be proved using the 
same approach as the one used in jT] to prove the corresponding statement for the Hahn 
ensemble. 

Let us fix a, 6 < 0 and 0 < /cq, go < 1- Assume that the conjecture is true, namely, that 
there exist such constants ci and C 2 that + C 2 N^^^u) converges to F 2 {u) as 

—> 00, where P 2 is the Tracy-Widom ditribution, q = q^ ■, k = koN, a = q^, and (d = q^. 
We write empathize the dependence on the parameters a, /?, g, N. Below we 

present a heuristic computation of the constants ci and C 2 , using a method suggested by 
Olshanski in [23| . 

Recall that the g-Hahn ensemble is the determinantal process with correlation kernel 
K{x, y) expressed through the monic orthogonal polynomials Pq = 1, Pi ,...: 

k-l 

K{x,y) = 2 Pi(a;)Pi(|/), x,y e X = {g"^|j = 0,..., - 1}, 

i=0 

where 

Pi{x) = '\/u{x)Pi{x). 

We will construct a selfadjoint operator D from the difference operator for the g-Hahn 
ensemble in such a way that: 

(i) 2) acts on T^(X); 

(ii) the correlation kernel K can be realized as a spectral projection P(A), here A is a 
certain part of the spectrum of the operator 2); 
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(iii) 2) formally converges to the Airy operator = g{u) — ug{u) as —>■ oo. 

In this approach the next step is to prove a stronger convergence of the operators and 
then deduce the convergence of the correlation kernels, see 0 , lisi- Unfortunately, we do 
not know how to do it rigorously at the moment. Nevertheless, in this way we can guess the 
scaling constants and check if they agree with numerical computations. 

Recall that the following g-difference equation holds j22| 

q-^{l - g")(l - al3q^^^)Pn{x) = B{x)Pn{x + 1) - (R(x) + D{x))Pn{x) + D{x)Pn{x - 1), 
where 

D{x) = a(3q{l - g^)(l - (5-^q^-^). 

The corresponding g-difference equation for Pn{x) is 

9-”(l - 9")(1 - aPq’'*')Pn{x) - B{x)A f + 1) 

V +1) 

- {B{x) + D{x))Pn{x) + D{x)j Pn{x - 1). 

Y 0J[X - 1) 

Let us denote by 2) the operator given by the right hand side. We are interested in 
the projection on the eigenvectors corresponding to n e {0,..., — 1}. Thus, the desired 

projection V is the spectral projection associated with the following segment: 

Since eventually we want to get the Airy kernel, which corresponds to the spectral pro¬ 
jection onto the positive part of the spectrum of the Airy operator, instead of T) we should 
consider 

2) = I) - g-(^-^)(l - g^"^)(l - a(3q^)3, 
where 3 is the identity operator. 

Let g{u) be a smooth function on M. Assign to it a function f{x) on j£ by setting f{x) = 
g{u), where x and u are related by x = ciN + C 2 N^^^u. 

Then we get 

f{x + 1) = g{u ± % g{u) ± 

We can hnd two values of ci (these corresponds to the rightmost and the leftmost particle) 

2 

and C 2 such that D ^ cN~i{^g{u) — ug{u)), where c is a constant. More precisely, 

_ log(p) 
log(go) 

and p is equal to 

{A^B^qo + A^Bqo + A^B + AB)K^ + {-2A^Bqo - 2ABqo - 2AB - 2A)K^ + {ABqo + Aqo + A+ 1)K ± 2qo\B5 
K{4:A^BqoK‘^ + {A'^B'^p^ - 2A^Bql - 2A^Bqo + _ 2A'^qo - 2ABqo + A^ - 2Aqo - 2A + 1)K + 4Aqo) ’ 
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where 


D = K{K^ -Kqo-K + qo){l - BK)A 

■ {A^B^K^qo - A^B'^K^qo - A^BK^qo - A^BK^ + ABKq^ + ABK + AK - 1) 

and A = q^, B = qQ, K = gg°. The constant C 2 can be expressed through a unique real 
solution of a cubic equation. We will not present here the formula for C 2 since it is quite 
lengthy. 

In the limit the eigenvalue equation D'ljj = A'0 turns into the equation = stfj after 

2 

the renormalization A = cN~3s. 

Next, we present some results of our computations. First, we plot the density function 

{Dk{ciN + C 2 uNs + 1 ) — Dk{ciN + C 2 uNi)) ■ C2N3 

for the same values of the parameters but different values of N. The blue graph is the density 
function of the Tracy-Widom distribution]^ 



Figure 3. Plot of the density function of the distribution of the rightmost 
particle with parameters N = 2000, q = 0.99, ko = 0.2, a = —1.1, b = —1.3, 
Cl = 0.84839, C 2 = 0.38999. 



Figure 4. Plot of the density function of the distribution of the rightmost 
particle with parameters N = 10000, q = 0.99, = 0.2, a = —1.1, b = —1.3, 

Cl = 0.84839, C 2 = 0.38999. 

*To plot the density of the Tracy-Widom distribution we used a table of its values http; //www. wisdom. 
weizmann.ac.il/~nadler/Wishart_Ratio_Trace/TW_ratio.html based on |5]. 
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We also plot E(TW) — E(qH^) in logarithmic coordinates, where TW is the Tracy-Widom 
distribution and qH^ is the distribution of the rightmost particle for the g-Hahn ensemble. 
The following pictures illustrate convergence of this graph to a line with slope — which 
agrees with our expectations m 



Figure 5. N = 10,20,... ,400; go = 0.99, = 0.3, a = -1.1, b = -1.3, 

Cl = 0.84839, C 2 = 0.38999. The blue line has slope — |. 



Figure 6. N = 10,20,..., 1000; go = 0.5, ko = 0.2, a = -1.1, b = -1.3, 

Cl = 0.69758, C2 = 0.47101. 
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